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Abstract 

We investigate the local integrability and linearizability of three di- 
mensional Lotka-Volterra equations at the origin. Necessary and suffi- 
cient conditions for both integrability and linearizability are obtained 
for (1,-1,1), (2,-1,1) and (1, — 2, l)-resonance. To prove sufficiency, 
we mainly use the method of Darboux with extensions for inverse Ja- 
cobi multipliers, and the linearizability of a node in two variables with 
power-series arguments in the third variable. 

1 Introduction 

In this paper we investigate the local integrability and linearizabihty of the 
origin for the three dimensional Lotka-Volterra systems, 

X = P = x{X + ax + hy + cz), 

y = Q = y{l^ + dx + ey + fz), (1) 
z = R = z{v + gx + hy + kz), 

where A, ^, 7^ 0. Lotka-Volterra systems are widely used with a diverse 
range of applications: population biology, chemical kinetics, laser physics, 
plasma physics, neural networks etc. [H 13 [6] 

We say that the system is integrable at the origin if there exists a change 
of coordinates 

X = x{l + 0{x,y,z)), y = y(l + 0(x,y,z)), Z = z{l + 0{x,y, z)), 

bringing ([1]) to a system orbitally equivalent to its linear system: 

X = XXm, Y = fiYm, Z = uZm, (2) 

where m = m(X, Y, Z) = 1 + 0{X, Y, Z). If the change of coordinates can 
be chosen so that m = 1 then we say the system is linearizable. 

We are interested in the case where the eigenvalues A, ^ and u are in 
the Siegel domain and have two independent resonances. Without loss of 
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generality, this means that we can assume (after a possible scaling of time) 
that X, fj^v £ 7, with gcd{X, fi, u) = 1, and that X, v > and fi < 0. We say 
that the origin has {X : fj, : z^)-resonance. 
Clearly, ([2]) has analytic first integrals 

(|) = X-^'Y^, and ^p = Y''Z-'', (3) 

which pull back to first integrals 

cp,=x-^y\l + 0{x,y,z)), and 02 = y"^"'^!! + 0(^, y, ^)), (4) 

of ([1]). Conversely, given two such first integrals, it is easy to construct a 
change of coordinates such that (j) and ip expressed in these new coordinates 
satisfy (I3D, and hence the transformed system is of the form ^ for some m. 

Restricting ([1]) to z = we can see the problem above as a generalization 
of the problem of classifying the integrability conditions for the system 

X = x{X + ax + by), y = y{n + dx + ey), X,fi£7,, Xfx<0. 

This problem was considered by several authors (Gravel and Thibault 
[18] , Christopher and Rousseau [10], Li et al. [19] and, for more general 
quadratic systems, by Zol§,dek [23], Fronville et al. [12] and Christopher et al. 
[9]) and gives a simple generalization of the Poincare center- focus problem. 
That is, to find conditions for which a local analytic first integral exists in a 
planar system. In this case where A + ^ = we have a complexified version 
of the classical center-focus problem. More recent work on integrability and 
linearizability of Lotka Volterra type systems with (p, — g)-resonance can be 
found in [Ifl [HI O [22] 

Recall that in the classical center-focus problem, there are currently only 
two known mechanisms for integrability: the existence of an algebraic sym- 
metry or the existence of a Darboux integrating factor. In the general case 
of p : —q resonance, however, other mechanisms appear to come into play: in 
particular, blow-down to a node and reduction to a Riccati equation [9[[T2]. 
For the Lotka- Volterra system, many of these conditions were subsumed in 
[1U| under a simple monodromy condition, applied to the neighborhoods of 
the invariant lines x = and y = 0, and the line at infinity. These conditions 
reduced essentially to finite checks on the nature of the singularities (finite 
and infinite) of the system. 

Our aim here is to see if a generalization of the Poincare center-focus 
problem to higher dimension gives a similarly simple list of integrability 
mechanisms. We give a complete classification of the integrability conditions 
for Lotka- Volterra equations with (1 : —1 : 1), (2 : —1 : 1) and (1 : —2 : 
1) resonant critical points at the origin. It was a surprise to us that the 
problems of integrability was much harder in this case, giving rise to new 
forms of argument which rely less on geometric properties than the form 
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of the power series concerned. This might be due to the fact that such 
resonant singularities mix the saddle and node-like properties of their two- 
dimensional counterparts. It would be interesting to know whether there 
were more geometric ways of obtaining the sufficiency of these conditions. 

Other work on 3D Lotka-Volterra equations has been done by Bobienski 
and Zolcidek [3j, who consider the finite singularity away from the axes 
planes, and give a number of mechanisms for the existence of a center in the 
{i : —i : A) case; Cairo and Llibre [H|5], who obtain a number of conditions 
for the existence of Darboux first integrals in terms of the parameters; and 
Basov and Romanovski [1], who take one of eigenvalues equal to zero. There 
has also been several works devoted to systems which are homogeneous 
(A = /U = J/ = 0) and hence reducible to a two-dimensional Lotka Volterra 
equation ([IH HZl EH [6]). 

We hope to extend this work to consider more general integrability phe- 
nomena in future work. In particular, the case where the system has a res- 
onant integrable critical point with non-integer ratio of eigenvalues should 
also be approachable using the methods here. In these cases, the first in- 
tegrals would no longer be analytic, but a product of powers of analytic 
functions: x°'y^z'^{l + 0{x, y, z)). 

The paper is organized as follows: In section 2, the Darboux method of 
integrability and inverse Jacobi multiplier together are explained as well as 
the relation between integrability and linearizability. In Section 3, we give 
the complete classification of integrability and linearizability conditions of 
system ([TJ with (1 : —1 : 1), (2 : —1 : 1) and (1 : —2 : l)resonance at the 
origin. 

2 Mechanisms for Integrability and Linearizability 

In this section we summarize the results we will need for understanding the 
integrability mechanisms which appear in our calculations. These methods 
can, in most cases, be extended to more general situations, but we focus 
here on the specific context of system ([T]). 

2.1 Reduction to the Poincare domain 

A singular point whose eigenvalues lie in the Poincare domain (that is, the 
convex hull of the eigenvalues does not contain the origin) can be brought 
to normal form via an analytic change of coordinates. In particular, a node 
with two analytic separatrices can have no resonant terms in its normal form 
and so must be analytically linearizable. 

We use this principle in two ways. Firstly, in many cases we can choose a 
coordinate system so that two of the variables decouple to give a linearizable 
node at the origin. If this is so, it just remains to find a linearizing transfor- 
mation for the third variable via some simple power series arguments. Sec- 
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ondly, and more rarely, we can perform a blow down to a three-dimensional 
system in the Poincare domain. Since this new system is linearizable, we 
can find two first integrals which we can pull back to first integrals of the 
original system. 



2.2 Darboux Integrability and Inverse Jacobi Multipliers 

Our second main mechanism for proving integrability is the existence of 
a Darboux first integral. The method is well know for two-dimensional 
systems [H [3 [HI [18j , but has been used also for higher dimensional systems 

Let 

ox oy oz 

be the associated vector field to the system ([T|). Given a polynomial F G 
C[x,y,z], a surface F = is called an invariant algebraic surface of the 
system ([T|), if the polynomial F satisfies the equation 

dF dF dF 
F = XF = P^ + Q^ + R^ = CfF (5) 
ox Oy oz 

for some polynomial Cj? G C. Such a polynomial is called the cof actor of 
the invariant algebraic curve -F = 0. One can note that from equation ([T]) 
that any cofactor has at most degree one since the polynomial vector field 
has degree two. 

To complete the study of integrals of parametric families, we will also 
need the notion of exponential factor which plays the same role of as an 
invariant algebraic surface in the case when two such surfaces coalesce. Let 
E{x,y,z) = exp{f{x,y,z)/g{x,y,z)) where f,g e C[x,y,z], then E in an 
exponential factor if 

XE = CeE, (6) 

for some polynomial C^; of degree at most one. The polynomial C^; is called 
the cofactor of E. 

A Darboux function is a function of the form, 

D = JJi^^'E^^o-^/^ 

where the Fi are invariant algebraic surfaces of the system, and E = exp{f /g) 
is an exponential factor. Given a Darboux function, D, we can compute 

Clearly, the function D is a non-trivial first integral of the system if and 
only if the cofactors Ci?. and Ce are linearly dependent. 
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For Darboux integrability in two dimensions, we seek a Darboux function 
which is either a first integral or integrating factor for the system. From the 
latter, it is possible to find a first integral by quadratures. 

In higher dimensions, the role of the integrating factor is taken by the 
Jacobi Multiplier. In the context of Darboux integrability, we usually con- 
sider the corresponding reciprocals: inverse integrating factors, and inverse 
Jacobi multipliers [2J. A function M is an inverse Jacobi multiplier for the 
vector field X if it satisfies the equation 

X(M) = Mdiv(X) div(X/M) = 0. 

A Darboux inverse Jacobi multiplier, D, must satisfy XiCp. + XqCe = 
div(X). 

In three dimensions, the existence of two independent first integrals im- 
plies the existence of an inverse Jacobi multiplier. Conversely, given just 
one first integral, (/>, and an inverse Jacobi multiplier, M, one can construct 
another first integral in the following manner. 

Suppose that the level surfaces (p = c are locally parameterized by some 
function z = fdx, y). Using the x and y coordinates to parameterize (p = c, 
we obtain a vector field 

d d 

P{x, y, fc{x, y))-Q^ + Q{x, y, fc{x, y))^- 

It can be shown [2] that 

defy 

M{x, y, fc{x, y)) — {x, y, /c(x, y)) 

is an inverse integrating factor for this vector field. Hence, by quadratures 
along (p = c, we can construct a second first integral ipc{x,y) for each value 
of c. The function ip(j,(x,y,z)ix,y) gives a second first integral of the system. 

However, in our case, we would like more control on the form of this 
second first integral. The following theorem allows us to give an explicit 
expression for the integral in certain cases. We use the usual multi-index 
notation X^ = x'^y^ to simplify the notation. 



Theorem 1. Suppose the analytic vector field 




has a first integral (p = x°'y^z'^{l + 0{x,y,z)) with at least one of a, /3, 
7 7^ and a Jacobi multiplier M = x^y^z^{l + 0{x,y,z)) and suppose that 
the cross product of {r — i — 1, s — j — l,t — k — l) and (a, (3,^) is bounded away 
from zero for any integers i,j,k > 0, then the system has a second analytic 
first integral of the form ip = x^'^y^''^ z^~^{l + 0{x,y, z)), and hence the 
system ^ is integrable. 
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Proof. Without loss of generality, we assume that a > 0. After an ana- 
lytic change of coordinates of the form {x,y,z) i— (x(l + 0(x, y, z)), y(l + 
0{x,y, z)), z{l + 0{x,y, z))), which will not alter the form of the vector 
field or the inverse Jacobi multiplier, we can assume that (f) = where 
6 = (a,/3,7). Furthermore, by absorbing the factor (1 + 0{x,y,z)) of M 
into the vector field itself, we can take the inverse Jacobi multiplier M 
to be , where 6 = {r,s,t). We take Aj = {Axi, Ayj, Azi) and write 

^(0,0,0) = (^^/^^i^)- 

From the hypothesis, we can take K > such that for all /, 

1(0-1-1) X 6\ > K. (7) 

Furthermore, since (/> is a first integral, then X(j) = gives for all / 

5-Ai = Q. (8) 

Since M is a Jacobi multiplier, we have X{M) = div(X)M, and writing 



where 1 = (1, 1, 1), we see that, for all /, 

{6 - I -!)■ Ai = (9) 

By hypothesis, {6 — 1 — 1) and 5 are linearly independent, and so ([8]) 
and ([9]) imply that there exists some constants kj such that 

Ai = ki{e-I-1) x6. (10) 

For ease of calculation we work with the associated 2-form Vt = P dy /\ 
dz + Q dz f\ dx + Rdx A dy rather than X. In this case, a function (j) \s a. 
first integral if and only if dc/) A = 0. Now 

Pdy Adz + Qdz Adx + Rdx Ady 



M XO 



dydz dzdx dxdy ^ vi-e+i 



y(kj(e-i-i) ^s).{^,^,^)x 

^ yz zx xy 



dx dy dz^s^ ^dx dy dz vi-d+i 

(11) 



ykA{e-i-i).(^, ^, ^)) A (5 . (^, ^, ^)) x^ 

^ X y z x y z 



ykj({e-i-i)-{—, ^, -) x'-'+^) A 

X y z 
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Thus, we have a formal first integral of the form 

I 

From d?]) and (fTO]l . we must have that \ki\ < K\Af\ for all I, and so ijj is in 
fact analytic. □ 



2.3 Integrability and Linearizability 

In some cases it is easy to deduce linearizability of a singularity from inte- 
grability. 

Theorem 2. // the system ([T|) is integrable and there exists a function 
^ = x'^y^z"'{l + 0{x,y,z)) such that X{^) = for some constant k = 
aX + + 7z^, then the system is linearizable. 

Proof To see this, note that (0) imphes that = ^X^^y^^Z^^ = 1 + 

0{x, y, z) satisfies Xi^l^^) = ^l^''{l-m). Hence (x, y, z) = {X^l ,¥^1 ,ZS,I) 
is a linearizing change of coordinates. □ 

The Lotka-Volterra equations have another property which was first 
noted in [10] for two dimensional systems. However, we do not use this 
condition explicitly in what follows. 

Theorem 3. Consider three dimensional Lotka- Volterra system ^ for which 
the three separatrices at the origin x = 0, y = and z = have cof actors L^, 
Ly and Lz respectively. If L^, Ly, and the divergence div{X) are linearly 
independent then the origin is integrable if and only if it is linearizable. 

Proof. Suppose that the origin is integrable. Then there exists two inde- 
pendent first integrals (p = x^f^y^ipi{x, y, z) and ip = y'^ z~'^ipi{x, y, z) where 
(x, y, z) = 1 + 0(x, y, z) and -01 {x, y,z) = 1 + 0{x, y, z) are analytic. The 
functions (pi{x,y,z) and ^i(x,y,z) obey the equations 

Xipi = ipiL^^, XV"! = V'l-^Vi' 

where 

Lipi = fJ-Lx ~ XLy, = —vLy + ^Lz- 

Since (p and ij) are first integrals, then dip A Q = 0, d^ A Q = and 
dip Adip = MQ where Q = P dy Adz + Q dz Adx + Rdx Ady as before, and 
-M is a Jacobi multiplier. One can easily show that 

where (/>(x, y, z) satisfies (/)(0, 0, 0) = — ^ / and (j) has a cofactor given by 
the divergence plus a linear combination of L^^ Ly and L^. Hence the co- 
factors La,, Ly, Lz and ^ are linearly independent. Note that the cofactors 
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L^^, and L^^ have no constant term. The condition on hnear indepen- 
dence imphes that we can find a change of coordinates X = x 93"^ 4'"^ > 
Y = y (ff^ ip^^ (f)^^ and Z = z tji^ (fp^ which hnearizes the system. □ 

3 Integrability and Linearizability Conditions 

In this section, we will give a complete classification of the integrability and 
linearizability conditions for the origin of ([1]) with (1 : — 1 : 1), (2 : —1 : 1) 
and (1 : — 2 : l)-resonance. We note that x = 0, y = and z = Q are always 
invariant algebraic surfaces with cof actors \ + ax + hy + cz^ ^ + dx + ey + fz 
and u + gx + hy + kz respectively. The necessary conditions were found by 
computing the conditions for the existence of two independent first integrals 
up to a given degree in the form (jH) using Maple. The degrees used 6, 10 and 
12 respectively. A factorized Grobner basis was then found using Reduce 
and finally the minAssGTZ algorithm in Singular [llj was used to check 
that the conditions found were irreducible. To prove sufficiency, we exhibit 
first integrals in the form 

For linearizability, we proceeded similarly: computing the conditions for 
the existence of a linearizing change of coordinates up to some finite order to 
find necessary conditions, and exhibiting a linearizing change of coordinates 
for sufficiency. In this case, the first integrals can be obtained easily by 
pulling back the first integrals of the linearized system 

3.1 (1 : — 1 : l)-resonance 

Theorem 4. Consider three dimensional Lotka-Volterra system ^ with 
(A,/x, i^) = (1,-1,1). The origin is integrable if and only if one of the 
following conditions are satisfied: 

1) ab — de = ac — 2ak + gk = ae + ah — de — eg = af + ak — dk — gk = 
bd + bg — de — dh = bf — ch — fh + hk = bk — ce + ek — hk = 
cd + eg — 2dk + fg — gk = ef — hk = 



2)6 = 


d = f = h 


= 




3)/ = 


g = h = b 


— e = 


d-a = 


3*) 6 = 


C = d=f 


-k = 


e-h = 


4)6 = 


c = d = f 


= k = 





4*)a = 


d = g = h 


= f = 


-- 


5) 6 = 


e = h = 
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Moreover, the system is linearizable if and only if either one of the conditions 
(2)- (5) or one of the following holds: 

1.1) a = c = d = f = g = k = 

1.2) a = bk — ch = d = e — h = f — k = g = 
1.2*) a-d = b-e = c = dh-eg = f = k = (} 

1.3) a — g = b — h = c— k = d — g = e — h = f — k = 

Proof. Cases 3*, 4* and 1.2* are dual to Cases 3, 4 and 1.2 under the trans- 
formation {x,y,z) !->■ {z,y,x), and do not need to be considered separately. 
The other cases are considered below. 

Case 1: If e 7^ 0, the system has an invariant algebraic plane i = 1 + ax — 
ey + kz = with cofactor L( = ax + ey + kz, and we have two independent 
first integrals 

(f>i = xy£~'^~'^, 4>2 = yz£~^~'^. 
When e = we have several cases. 

i) b = h = 0: This is a subcase of Case 5. 

ii) h = 0, b 7^ 0: In this case, we have a = f = k = g + d = and we get 
an exponential factor £ = exp(dx — by + cz) with cofactor dx + hy + cz, and 
first integrals (f)i = xy£'~^ and (j)2 = yz. 

iii) b = 0, h 7^ 0: In this case, we have a = d = k = and c = f{b/h — 1). 
We get an exponential factor £ = exp{gx—hy+fz) with cofactor gx+hy+fz. 
This gives first integrals (f)i = xy and (^2 = yz£~^. 

iv) b, h 7^ 0: In this case, we have a = k = 0, c = f{b/h—l) and g = d{h/b— 
1). We get an exponential factor £ = ex.p{dhx — bhy + bfz) with cofactor 
dhx + bhy + bfz. This gives first integrals (f)i = xy£~h and 02 = yz£~b . 

Case 1.1: If e 7^ 0, we have an invariant plane 1 — ej/ = 0, and the change 

of coordinates {X,Y,Z) = {x{l — ey)~^ ,y{l — ey)~^ , z{l — ey)~^) linearizes 
the systems. When e = 0, we replace (1 — ey)~e and (1 — ey)~e above by 
exp{by) and exp{hy) respectively. 

Case 1.2: When e 7^ The linearizing change of variables is given by 
{X,Y,Z) = {x{l - ey + fz)--e ,y{l - ey + fz)-\zil - ey + fz)-'). When 
e = 0, then either 6 = and / 7^ 0, and we have a linearizing change of 
coordinates {X,Y,Z) = (.t(1 + fzY'^,y{l + fz)-^,z{l + fz)-'^), or / = 
and we linearize by {X, Y, Z) = (xexp(6y — cz),y, z). 

Case 1.3: In this case we have an invariant plane £ = 1 + ax — by + fz 
with cofactor Li = ax + by + fz and the linearizing change is {X, Y, Z) = 
{x£-^,y£-^,z£-^). 
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Case 2: In this case we have an invariant algebraic plane i = 1 — ey = with 
cofactor Li = ey. Writing Y = y/(l — ey), we obtain Y = —Y. Furthermore, 
the equations 

X = x{l + ax + cz), z = z{l + gx + kz), (12) 

gives a linearizable node. Thus, there exists X = x{l + 0{x, z)) and Z = 
z{l + 0{x, z)) such that X = X and Z = Z. 

Case 3: The transformation X = x/{l + ax — by), Y = y/{l + ax — by) 
gives 

X = X{l + cz- acXz), Y = Y{-1 - acXz), z = z{l + kz). 

The first and third equations give a linearizable node, and so we can find 
a change of coordinates X = X{1 + 0{X,z)) and Z = z{l + 0{X,z)) 

such that X = X and Z = ^ (in fact Z = z/{l + kz)). Since Y = 
y(— 1 — acX{X , Z) z{Z)), it is sufficient to find a function £{X,Y) such 
that i{X,Z) = X{X,Z)z{Z), then the substitution Y = Ye""'' will give 

Y = — y, and the system is linearized. 

Writing X{X, Z)z{Z) = Ei+j>o aijX^'Z^, it is easy to see that 1{X, Z) = 
Si+j>o i^^ij-^^^'' gives a convergent expression for i. 

Case 4: This system has two invariant algebraic planes 1 + ax = and 
1 — ey = with cofactors ax and ey respectively, and a linearizing change 
of coordinates 

(X, Y, Z) = (x(l + ax)-\y{l - ey)-\z{l + ax)"^ (1 - ey)"^). 

When a = 0, we have an exponential factor exp(a:), and we replace 
(1 + ax)^/" by ex.p{gx). Similarly, when e = 0, we have an exponential 
factor exp{y) and replace (1 — ey)^l^ by exp(— /iy). 

Case 5: The equations in x and z are independent of y, and give a lineariz- 
able node 

X = x(l + ax + cz), i = z{\ ^ gx ^ kz). 

Hence there exists a change of coordinates X = x(l + 0{x,z)) and Z = 
z{l + 0{x, z)) such that X = X and Z = Z. The remaining equation is 
given by 

y = y{-l + dx{X,Z) + fz{X,Z)). 

Suppose there exists a function i{X, Z) such that ({X, Z) = {dx{X, Z) + 
fz{X, Z)), then the transformation Y = ye~^ will gives Y = —Y. 

Writing dx{X, Z)+fz{X, Z) = Ei+,>o aijX'Z^, we have £ = Ei+,>o 
which is clearly convergent. 

□ 
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3.2 (2 : — 1 : l)-resonance 

Theorem 5. We Consider three dimensional Lotka-Volterra system ^ with 
(A,/x,i^) = (2,-1,1). Under the following (l)-(ll) conditions, the origin is 
integrable: 

1) ef — hk = ah — ah — de + eg = ac — 3ak + 2gk = ae + ah — de — eg = 
af + ak — dk — gk = hd + bg — de — 2dh + eg = bf — ch — 2fh + 2hk = 
bk-ce + 2ek - 2hk = cd + eg - 3dk + 2fg - gk = 

2) b + e = c = d = f = k = 

3) b + h = c + k = d = e-h = f-k = 

4) a-d = b-e = f = g = h = 

5) b + e = d = f = h = 
Q)a = d = f = g = h = 

7) b = c = e-h = f-k = 

8) b = e = h = 

9) b = f = h = 
I0)b = c = f = k = 

n)b = c-Ak = e + h = f + k = 

Moreover, the origin is linearizable if and only if the system satisfied either 
one of the conditions (2)-(10) or one of the following conditions: 

1.1) a = c = d = f = g = k = 

1.2) a — d = b — e = c = dh — eg = f = k = 

1.3) a = bk — ch = d = e — h = f — k = g = 

1.4) a — g = b — h = c — k = d — g = e — h = f — k = 

Proof. We treat the cases individually. 

Case 1: If e 7^ 0, then we have an invariant surface i = 1 + |x — ey + /cz = 
with cofactor Li = ax+ey+kz. From this we can construct two independent 
first integrals 

Q b_ -1 h 

= X y i e and 4'2 = y e . 

When e = and /i 7^ then we have an exponential factor d. = exp((d + 
g)x/2 — hy + fz) with cofactor [d + g)x + hy -\- fz with corresponding first 
integrals 

Q b_ 1 

(pi = xy I h and (p2 = y zi . 

If e = and h = then we can assume 6 7^ since otherwise we are in 
Case 8. The conditions then give an exponential factor £ = exp{gx + by — cz) 
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with cofactor 2gx — by — cz and corresponding first integrals ipi = xy"^ £ and 
(j)2 = yz. 

Case 1.1: If e 7^ 0, the Hnearizing change of coordinates is given by 

{X,Y,Z) = {xil - ey)--^,y{l - ey)-\ z{l - ey)-^). 

When e = 0, we have an exponential factor exp(y), then replace (1 — ey)~e 
and (1 — ey)^e by exp(6y) and ex.p{hy) respectively. 

Case 1.2: For 6 7^ 0, a linearizing change of variables is given by 

(X, Y, Z) = + |x - hy)-\y{l + - by)-\z{l + ^x - by)-^). 

When 5 = 0, either h = and a 7^ 0, and we can linearize by 

{X,Y,Z) = (x(l + ^z)-\y{l + ^z)-\z{l + \z)-^l^\ 

or a = and we linearize by (X, y, Z) = (x, y, z exp(6?/ — gx)). 

Case 1.3: This case is exactly similar to case 1.2. 

Case 1.4: In this case we have an invariant surface i = 1 + ^x — by + fz with 
cofactor Li = ax + by + fz and the linearizing change is give by {X, Y, Z) = 

{xe-^,ye-^,z£-^). 

Case 2: In this case we have invariant planes ii = 1 + ^x + by = 0, 
£2 = 1 + — ^xy = and £3 = 1 + by = with cofactors L^^ = ax — by, 
Li^ = ax and Li^ = —by respectively. 

— - — 

The change of coordinates {X,Y,Z) = {x£^ ^hv^s '-^^2 "^3 ) linearizes 
the system. When 5 = we can replace £3''^ by exp(/iy) and, when a = 0, 
we replace £2^^"" with exp(— x(l — by)g/2). 

Case 3: After the change of coordinates {X, Y, Z) = {x,xy, xz), the system 
becomes 

X = 2X + aX^-eY-fZ, Y = Y{l + aX), Z = Z{3 + {a + g)X). (13) 

The critical point at the origin of ()13p is in the Poincare domain and hence 
is linearizable via an analytic change of coordinates which can be chosen to 
be of the form {X,Y,Z) = {X - eY + fZ + 0{2),Y{1 + 0{1)), Z{1 + 0{1))). 

The two first integrals (p = X'^Y"^ and ■i/' = X~^Y Z of the linear 
system pull back to first integrals of the form 

=xy2(i + 0(i)), and 4^2 = yz{l + 0{l)). 



12 



The expression ^ = x°'y°'~^^z ^ satisfies ^ = {a — 2g)^ and so the system is 
hnearizable from Theorem [2] 

Case 4: The transformation X = x / {1 + ax /2 — by), decouples the terms 
in X and z to give 

X = X{2 + {1- aX/2)cz), z = z{l + kz), (14) 

with a hnearizable node at the origin. Thus, there is a linearizing change 
of coordinates X = X{X, z) = x{l + 0(1)), Z = Z{z) = z{l + 0(1)) which 

brings (HH) to the form X = 2X and Z = Z. 

Taking Y = y/{l + ax/2 - by), we get Y = Y{-1 - acXz/2). It is suf- 
ficient, therefore, to find e{X,Z) such that e{X,Z) = X{X,Z)z{Z). Then 
the substitution Y = Ye"'^^'^ will linearize the system. 

lix{X,Z)z{Z) = Y.^+,>o'^^,X'Z'J, then i{X,Z) = J^^^^^^^X^Z^ 
gives a convergent expression for i. 

Case 5: A change of coordinates {X,Y,Z) = {j^,y,i:^), brings the 
system to the form 

X = X{2+aX+cZ){l+bY), Y = -Y{l+bY), Z = Z{l+gX+kZ){l+bY). 

(15) 

After rescaling the system by (1 + bY), we have Y = —Y and the first and 
third equation give a hnearizable node. This implies the original system is 
integrable. However, since Y = —Y the system must also be hnearizable by 
Theorem [2j 

Case 6: A linearizing change of coordinates is given by 

(X, Y, Z) = (x(l - ey)-|(l + kz)~i,y{l - ey)-\z{l + kz)-^). 

When A; = 0, we have an exponential factor exp(z), and we replace {l+kz)~k 
by exp(— cz). Similarly, when e = 0, we have an exponential factor exp(?/) 
and replace (1 — ey)~e by exp(6y). 

Case 7: The system has an invariant algebraic plane i = l + ax/2 = with 

-I -1 d — g-\-a 

cofactor = ax yielding a first integral (j) = x" y^ z i ^ . We also have 
an inverse Jacobi multiplier IJM = x'2y^i^2^~^i . 

d—g+a 

When a = 0, an exponential factor exp(x) will appear and replace I °- 
and ^~2~"S"+f by exp(^^x) and exp((| — d)x) respectively. Theorem 1 
therefore guarantees the existence of a second first integral of the form 
^ = x~'^l'^y~'^z{\ + 0(1)). From these two integrals it is easy to con- 
struct integrals of the form required as <f)\ = (lP"il)~'^ = xy^(l -|- . . .) and 
<^2 = 4'^i/j~'^ = yz{l -|- . . .). Since ^ = xi~^ satisfies ^ = 2^, the system is 
linear by Theorem [2l 
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Case 8: In this case, our system is 

X = x{2 + ax + cz), y = y{—l + dx + fz), z = z{l + gx + kz). 

The first and third equations gives a hnearizable node, so it suffices to find 
e{X, Y) such that 

i{X,Z) = dx{X,Z) + fz{X,Z) (16) 

The transformation Y = ye~^ then hnearizes the second equation. 

Writing dx{X,Z) + fz{X,Z) = T.i+j>oaijX' Z^ , we obtain e{X,Z) = 
Yli+j>o ^+7^*-^''' which is clearly convergent. 

Case 9: The system ([1]) can be written as 

x = x{2 + ax + cz), y = y{—l + dx + ey), z = z{l + gx + kz). (17) 

The first and third equations in (|17|) give a linearizable node. To linearize 
the second equation, we seek an invariant surface of the form H. + xn = ^ 
with cofactor dx + ey where i = £{X, Z), x = xi^j ■^)i ^(0) = 1. The 
change of variable Y = will then linearize the second equation. 
To find £ and x we therefore need to solve 

X-x = e£, e = dx£. (18) 

To find £, we write £ = and solve ip = dx. Let ip = ^j_|_j>o '^ij^^^'' ^ 
then 

{2i + j)cijX'Z^ = dx{X, Z) = dX+ J2 dijX'Z^, 
i+j>0 «+i>i 

for some dij. Clearly, cio = cqi = 0, Cij = for i + j > 1. The 

convergence of X^j_|_j>i dijX^Z^ , guarantees the convergence of -0 and hence 
£. Furthermore, it is clear that £ will contain no term in Z. 

Now, writing £ = Y^bijX^Z^ , and noting that oqi = 0, we find that 
X = 2i+j-i '^V'^^ ■^^ gives a convergent expression for x- 

Case 10: In this case the system ([TJ reduces to 

X = x{2 + ax), ij 

The transformation 

X 

would linearize the first and the second equations in ()19p if £ + = were 
the defining equation of an invariant algebraic surface with cofactor dx + ey 
where £ = £{X), x = x{^) and £{fi) = 1. This condition is equivalent to 

X-X = e£, i = dxi. (21) 



y{—l + dx + ey), z = z{l + gx + hy). (19) 



l + ax/2' 



Y 



+ xy 



(20) 
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The second of these equations is clearly satisfied if we take I = = 
(1 - aX/2)-2/'^ (when a = let £ = = e^). Writing x = E«i^* and 
i = Yl biX^ ^ equation ([2T]) is satisfied if we set Oj = bi/{2i — l). The resulting 
function x is clearly convergent. 

To linearize the third equation, it is suffices to find "y{X, Y) such that 
7 = g x{X) + hy{X,Y), and take Z = zexp{—j). 

Writing, 

7= Yl c.,^'^^ x{X) = J2diX\ y{X,Y)= J2 ^ijX'Y'^ 

i+j>0 i>0 i+j>0 

we see that we require 

Y,{2i - j)ci,X'Y^ = Y,9 d^X' + Y,h Ci.X'Y^ = ^ h.X'YK 

If y{X,Y) contains no terms of the form X^Y"^^ , then we can set Cij = 
{gdi + heij)/{2i — j) for 2i ^ j, and Cij = otherwise, to find a convergent 
expression for 7. 

Therefore, it just remains to show that the inverse transformation y = 
y{X,Y) from ^ contains no term like (XY^)"-. Prom 

suppose n + 1 = 2m for some m. It is suffices to show that £x^ contains no 
term X"^ or in another word ^x^™"^ has no term X™. Since 

X = ^ = 2X^ 
^ dt dX' 

then, from equation ()2ip . we obtain 

However, the term in X™ in x?"^ clearly vanishes on the left hand side of 
()22p . so that either e = 0, in which case x = 0, or the coefficient of X^ in 
T^2m-i vanishes. Thus y in (|2Up contains no (Xy^)" and we have established 
the existence of a linearizing change of coordinates. 

Case 11: The transformation w = yz brings the system to the form 

X = x{2 + ax — 4/z), w = wxd, z = z{l + gx — fz) — ew, (23) 

where d = d + g. In this case we seek and expression = w^ipiix, z) such 
that ip = X. If such a tp exists, then (j) = we~'^'^ is a first integral. 
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We write the vector field as X = Xq + + X^ , where 

d d d ~ d 

Xo = x{2+ax-4:fz)- — \-z{l+gx-f z) — , Xi = -ew—, X^ = wxd—. 

ox oz oz ow 

Prom X^ = X, we get: 

Xo'il^o = x, Xoil)k + kxdil:k = -Xiil)k-i {k>{)). (24) 

We now solve the equation ()24p . To do this, we first show that for 
every B = Yli+jyo^ij^^ ^ there exists an A = Ylii+j>o^ij^^'^'^ ■> such that 
(Xo + kdx) A = B. Since 

{XQ+kdx)A = ^ {2i+j)aijx'z^+ ^ {ia+jg+kd)aijx'^'^ z^ - ^ {4i+j)faijx''z^ 

i+j>0 i+j>0 i+j>0 

we find that the aij must satisfy 

{2i + j)aij + ((i - l)a + jg + kd)ai-ij - (4i +j- l)/aij_i = kj. 

Now in (j24p we can solve the equations term by term provided that 
the right hand side of the equations have no constant term. However, this 
follows from the stronger observation that we can choose ipi in ()24p to be 
divisible by x. To show this, we can suppose by induction that x divides 
the right hand side of (j24p (for A; = this is immediate). This implies that 
z{l+gx — fz)^^ is divisible by x and this in turn shows that is divisible 
by X. Writing ipkix, z) = g{z) + x h{x, z), we see that g'{z) = 0, so that g is 
a constant. Clearly ipk — 9 also satisfies and we proceed by induction. 

There is therefore no obstruction to solving these equations ()24p recur- 
sively, and standard majorization techniques imply that the resulting series 
will be convergent. 

Thus, we find a first integral (pi = yze~'^^ of the original system. The 
system also has an inverse Jacobi multiplier 

3 1 I 9 + °/^ g + a/2 
X2y d Z d . 

Theorem 1 therefore guarantees a second first integral 

1 g + a/2 , 9 + a/2 
'll} = X 2y ~d-z ^~(l + 0(l)), 

from which we deduce the following first integral 

2d— a 

(p2 = (pi' = xy\l + 0{l)). 
When d = then x^l'^y is an IJM and we proceed as before. 

□ 
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3.3 (1 : — 2 : l)-resonance 

Theorem 6. The origin of the three dimensional Lotka-Volterra system m) 
with {\,^,u) = (1,-2,1) is integrable if and only if on of the following 
conditions holds: 

1) ab + ah — de — eg = ac — 2ak + gk = 

ae + 2ah — de — 2eg = af + 2ak — dk — 2gk = 

bd + 2bg — de — dh — eg = bf + ce — 2ch — ek — fh + 2hk = 





bk — ce 


1 7 7 7 7 J C\ / "\ 7 7 /* 

+ ek- hk = cd + 2cg - 2dk + fg - 


2gk 




ef + ek 


- 2hk = 




2) 


a = d = 


f+k=g=h=0 




3) 


a-2g-- 


7 1 O J T C I 1 1 

= b + e = c — 3k = a = f + k = n = 





4) 


c = d = 


f = 9 = 




5) 


a = d = 


f = 9 = 




5*) 


c = d = 


f = k = 




6) 


b-h = 


d = f = 




7) 
8) 


a- g = 
b + h = 


c-k = d = f = 

c + k = d = e — h = f — k = g = 




9) 


b + 3h = 


= c-3k = d = e + 2h = f + 2k = g 


= 


10) 


b = c = 


d = e-h = f -k = 




10*) 


a — d = 


b- e = f = g = h = 




11) 


b = e = 


h = 




12) 


a + d = 


b = cd-2dk-gk = f + k = h = 




13) 


a — d = 


b-e=f+k=g=h=0 




13*) 


a + d = 


b = c = e — h = f — k = 




14) 


c-2k = 


= d = f + k = g = h = 




14*) 


2a- g-- 


= b = c = 2d + g = f = 




15) 


a + g = 


b + 3h = c-3k = d = e + 2h = f^ 


- 2k 


16) 


3a- g ^ 


= b = 2c-k = 3d + g = e + h = f = 


= 


17) 


a + g = 


b+h=c=d+g=e+h=f=0 




18) 


3a- g - 


= 3b + h = c = 3d + 2g = 3e-2h = 


/ = 


19) 


3a- g-- 


= 3b + h = c + k = 3d + 2g = 3e-2h = 


20) 


b = d = 


f+k=g=h=0 




20*) 


a + d = 


b = c = f = h = 




21) 


a + d = 


b = c = f = k = 
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Moreover, the system, is linearizahle if and only if either one of the conditions 
(2)-(8), (lO)-(U), (16)-(n), (20)-(21) or one of the following holds: 

1.1) a = bk — ch = d = e — h = f — k = g = 
1.1*) a — d = b — e = c = dh — eg = f = k = 

1.2) a — g = b — h = c— k = d — g = e — h = f — k = 

Proof Cases 1*, 5*, 10*, 13*, 14* and 20.1* are dual to Cases 1, 5, 10, 13, 
14 and 20.1 under the transformation (x, z) i-> {z, x), and we do not consider 
them further. The other cases are considered below. 

Case 1: If e 7^ 0, the system has an invariant algebraic plane i = 1 + ax — 
+ kz = with cofactor Li = ax + ey + kz and produces two independent 

first integrals = x"^ y and = y z^ . 

If e = 0, we have several sub cases: 

i) k 7^ 0: We have 6 = ^ = and the system has first integrals 01 = 
x'^y£~^^^^ and (t>2 = y z'^ i'''^~k . 

ii) h 7^ 0: We have a = k = and we have an exponential factor i = 
exp{ — {d + 2g)x + hy — fz) with cofactor = —{d + 2g)x — 2hy — fz which 

yields first integrals (pi = x"^ y £h and 4>2 = y 

iii) h = k = 0: If 6 = we are in Case 11, so we assume that 6 7^ which 
implies that a = d + 2g = f = {). Then there exists an exponential factor 
a. = exp(2gx + by — 2cz) with cofactor = 2gx — 2by — 2cz which yields 
two first integrals 0i = x'^yi and (p2 = yz^- 

Case 1.1: If e 7^ 0, the change of coordinates {X,Y,Z) = (x(l — |y + 

fz)~~^,y{l — |y + fz)~^,z{l — |y + fz)^^) linearizes the system. When 
e = 0, then either fc = or 6 = 0. In the first case, we can linearize via 
{X,Y.Z) = {x exp(^6y — cz),y,z), and in the second, taking 7^ 0, we can 
linearize via {X,Y,Z) = {x{l + kz)~^ ,y{l + kz^^ , z{l + kzy^). 

Case 1.2: In this case we have an invariant plane i = 1 + ax — ^y + cz 
with cofactor = ax + by + cz and the linearizing change is {X, Y, Z) = 



Case 2: The system has an invariant algebraic plane ^ = 1 — |y — /2; = 
with cofactor = ey — fz. The substitution 



{xe-^,ye-^,ze-^). 



Y = yr'{l-fzf, 



z = 



z 



(25) 




xe 



will 



(26) 



18 



An obstruction to the existence of (p is possible only \i y = y{Y, Z) in 
(j26p contains a term of the form {Y Z'^)'"'. According to (j25p . one can find 

^-i + fy(i + /z)2-4^/-2) ^(1 + /^) ' 

and hence y contains no term of the form (yZ^)". 

Case 3: If 51 = then we obtain a sub-case of Case 2 and if / = we 
obtain a subcase of Case 5*. Hence, we shall assume that / and g are both 
non-zero. The system has two first integrals: 

^ xyjgexy - 2gx + 2fz - 2fz'^) 
^ {l-2gexy + 2gx-2fz + pz^y 

and 



(1 - 2gexy + 2gx - 2fz + pz'^ f ' 
From these, we obtain two first integrals of the desirable form 



and 



H = -;^ + -VV = a^'y(i + ---) 



-2 



The expression ^ = xyz satisfies C = and so the system is lin- 
earizable by Theorem [2j 

Case 4: The system ([T]) reduces to 

X = x{l + ax + by), y = y{—2 + ey), z = z{l + hy + kz). (27) 

The change of coordinates Y = y/{l — ey/2) gives Y = —2Y. To linearize 
the first equation, we seek an invariant algebraic surface A(Y) + B(Y)x = 
with cofactor ax + by, where A and B are analytic with A{0) = 1. Thus we 
seek A and B to satisfy the equation 

-2Y A'{Y) = byA{Y), -2Y B'{Y) + B{Y) = aA{Y). (28) 

The first equation gives ^ = (1 - ey/2)''/^ (or exp(5y/2) when e = 0). 

Writing A = "^iXjCLiY"^, we have B = "^i^QCLj^^Y^, which is clearly 
convergent. 

The substitution X = x/{A + Bx) linearizes the first equation of ()27p . 
In the same way, we can find an invariant surface C(Y) + D{Y)z = with 
cofactor hy + kz to so that Z = z/{C + Dz) linearizes the third equation of 
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Case 5: In this case the system ([T]) becomes 

X = x{l +by + cz), y = y{-2 + ey), z = z{l + hy + kz), 

The substitution Y = y/{l — ey/2) hnearizes the second equation. As in 
the previous case, we can find an invariant surface C(Y) + D{Y)z = 
with cofactor hy + kz where C and D are analytic with C(0) = 1. Then 
Z = z/{C + Dz) Hnearizes the third equation. Such a C and D must satisfy 

-2YC'{Y) = hyC{Y), -2Y D' {Y) + D{Y) = kC{Y). (29) 

We seek a transformation X = xe~^ which linearizes the first equation. 
Such a (j) satisfies 

4>{Y,Z) = by{Y) + cz{Y,Z). (30) 

Taking (p = ^ aijY^Z^ , we find that (f> = ^{2i — j)aij. It is clear that (j30p 
can be solved analytically for (p as long as z(Y, Z) contains no term of the 
form {YZ'^Y. Since 



)(y) 

^ ' i>0 

it is sufficient to show that CD^"^^ contains no term in y". 
However, from the defining equation ()29p . we have 

CD-n-. = U-2Y^^ + D) D^-^ = U--^-^ + D^- 



k \ dY J k\n dY 

which clearly does not contain any term of degree y". When A; = 0, we have 
D = from (j29p and the corresponding result is trivial. 

Case 6: In this case, the change of variables X = xm(y) and Z = zm{y) 
for some analytic function m{y) with m(0) = 1, brings the system to the 
form 

X = X{l + ax + cz + hy + —y{-2 + ey)), 

m 

y = y{-2 + ey), (31) 

Z = Z{1 + gx + kz + by + —y{-2 + ey)). 

m 



We choose m so that 



1 + by + —y{-2 + ey) = -. 

m m 



This equation has an explicit solution, 

m = - e?//2)-(e+i) J y~i(l _ ey/2)t-i 
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which can be seen to be analytic in y with m(0) = 1, taking the definite 
integral as a Laurent series in odd powers of ^yy. After scaling by m{y), the 
system becomes, 

X = X{l + aX + cZ), y = y{-2 + ey)m{y), Z = Z{1 + gX + kZ). 

Clearly the first and third equations gives a linearizable node. The second 
equation can be linearized by a substitution Y = i{y), such that 

y(2-ey)m{y)^ = 2e{y), £{0) = 0, £'(0) = 1. 

This is clearly solvable analytically once we have determined m(y) as above. 

Case 7: In this case the system has an invariant algebraic plane i = 1 — 
ey/2 = with cofactor = ey producing a first integral 

(j) = X Z i e 

It is easy to see that we have an inverse Jacobi multiplier 

IJM = - ey/2)e-^+^ (32) 

When e = 0, we take (p = x exp(^y) and IJM = yi exp(^^y) 

above. By Theorem 1, we can obtain a second first integral tp = xy~^z~'^ (l+ 
0(1)). Now we can construct two independent first integrals of the desired 
form as 

(pi = (p'^ip"'^ = x'^y{l + m{x,y,z)) and (p2 = (p'^ip~'^ = yz'^{l + m{x,y, z)). 
In this case the linearizing change is given by 

{X, Y, Z) = {x{l + m)^ (-2 + ey)^,y (-2 + ey)-^,z (1 + m)^ (-2 + ey) 5 ). 



Case 8: This case has invariant algebraic surfaces ii = 1 + ^y — cz = and 

£2 = l + ax + abxy — ^xz = 0, with cofactors L^^ = —by — cz and L^^ = ax, 
which allow us to find two independent first integrals: 

(pi = x'^ y £±£2'^ and (p2 = yz^£'^^, 

and linearizing change of coordinates {X,Y,Z) = {x £±£2^ ,y £1^ , z £^^). 

Case 9: In this case £ = ax{l + ^z) — ahxy + {1 + hy + kz)"^ = is an 
invariant algebraic surface with cofactor = ax — 4hy + 2kz, giving rise to 
the first integrals 

(pi=x'^y£~'^ and (p2 = yz'^. 
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Case 10: The change of variables {Y, Z) = ( ^^j^^g^/g ' i+fz-ey/2 ) S^-'ves a 
new system 

x = xil + ax), Y = -2Y{1 + fgxZ), Z = Z{1 + gx - fgxZ), 

The first and the third equations obviously gives a node and therefore there 
is a linearizing change of coordinates X = X{x), Z = Z{x, Z). To linearize 
the second equation, it is suffices to find ip{X, Z) such that ^ = fgxZ 
and use the transformation Y = Ye^^. Setting ip{X,Z) = Y^^^-^^aijX^ Z^ 

and x{X , Z)Z{X ^ Z) = "^bijX^Z^ , we find that Uij = dij/{i + j), giving a 
convergent expression for -0. 

Case 11: Then the system reduces to 

X = x(l + ax + cz), y = y{—2 + dx + fz), z = z{l + gx + kz), 

The first and third equations gives a linearizable node. We denote the 
linearizing coordinates by X and Z. 

The transformation Y = ye~'^ will linearize the second equation if cf) 
can be chosen so that (p{X,Z) = dx{X,Z) + fz{X,Z). Let (f){X,Z) = 
^i+j>o(^ij^'^^ ™d dx{X,Z) + fz{X,Z) = Y^i+j^ohjX'-Z^ , then we re- 
quire Uij = hij/{i + j), which gives a convergent expression for (p. 

Case 12: The system (P) reduces to 

X = x{l — dx + cz), y = y{—2 + dx + ey — kz), z = z{l + gx + kz). 

The system has invariant algebraic surfaces ii = 1 — dx — ^y + kz = and 
£2 = 1 — ^y — ^xy + ^yz = with cofactors L^^ = —dx + ey + kz and 
= ey. The substitution Y = y£i £2'^ linearizes the second equation, and 
the first and third equations define a linearizable node. 

Case 13: In this case system ([T|) becomes 

X = x{l + ax + by + cz), y = y{—2 + ax + by — kz), z = z{l + kz), 

The system has an invariant algebraic plane £ = 1 + kz = with cofactor 
Li = kz. By mean of a change of variables 

iX,Y,Z) = {x{l + kz)-^,y{l + kz),z{l + kz)-'^) 

we arrive to the system 

X =X{1 + aX{l - kZ)^T^ + bY{l - kZ)), 

Y =Y{-2 + aX{l - kZ)-i + bY{l - kZ)), (33) 

Z =Z, 
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with first integral (j) = X ^YZ^ and inverse Jacobi multiplier IJM = X^Y. 

We cannot apply Theorem 1, as there are non-negative integer values of 
i, j and k for which the cross product of (1 — i, —j, —1 — k) and (—1, 1, 3) 
is zero. However, this is only true when, (1 — i, —j, —1 — k) = a (—1, 1, 3) 
for some a. Clearly, the only possibility is when a = — l,i = 0, j = l and 
k = 2. 

However, in this case, the proof of Theorem 1 will still work as long 
as ^(0,1,2) = because, in this case, (fTO]l will still hold. But it is clear 
that in (|33p that there are no terms in YZ^ in the cofactors of X, y or 
Z, so indeed ^(0,1,2) = ^-nd we have a second first integral of the form 
ijj = X~^Z{1 + 0{X,Y, Z)). We get first integrals in the required form by 
pulhng back the first integrals (pi = (l)/'tp^ = X'^Y{1 + 0{X, Y, Z)) and (t)2 = 
(plijj = YZ'^{1 + 0{X,Y, Z)) to the original coordinates. The substitution 
{X,Y,Z) = {Z / ijj , (1)2 / Z"^ , Z) linearizes the system. 

When A; = 0, we replace (1 + kz)~k by exp(— cz) and proceed similarly. 

Case 14: In this case we have two invariant algebraic planes ^1 = 1 — |7/ — 
fz = and £2 = ! — fz = with cofactors = ey — fz and L^^ = —fz. 
This allow us to construct a first integral (p = yz^^i^ ^-nd inverse Jacobi 

2—- -—1 

multiplier IJM = £^ " ^1 • Now, Theorem 1 guarantees a second 

first integral ijj = x~^yz'^{l + 0{x, y, z)). We write the first integrals in the 
desired form as 

(j)i = = yiyl+m^x^y, z)) and (t)2 = (j) = y z^{l+m{x,y, z)) . 

The following change of coordinates linearizes the system 

(X,y,Z) = (x(l + m)l(l-/z)-i(l + m)-^y(l-/z)2(l + rn),z(l-/z)-i). 

Case 15: In this case li = {1 + hy + kz^ — kgxz = and ^2 = gx{—2 + 
gx + 2hy — 2kz) + (1 + /ly + kz)'^ = are invariant algebraic surfaces with 
cofactors L^^ = —Ahy + 2kz and L^^ = —2gx — Ahy + 2kz. The two first 
integrals are given by 

(Pi=x'^yq^q^ and (p2 = V q 

Case 16: The system has invariant algebraic surfaces ii = 1 + 2ax + 2cz + 
a^x^ — 2ceyz = and £2 = 2a?'x^ + 2cz — ceyz + 2ax = with cofactors = 

_3 

2ax + 2cz and L^^ = 1 + 2ax + 2cz. This gives a first integrals (j) = xy zi-^ ^ 
and ijj = x~^i2£i ^^"^ , though the latter integral is not of the required form. 
However, when z = 0, we have '0 = 2a. Hence, = {ip — 2a)x/{2cz) = 1 + 
0(1) is analytic and satisfies ^ = ^(— 3ax+ey — cz). From this we construct a 
linearizing change of coordinates {X,Y,Z) = {xi^ ^^"^ ,yl^ 
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When c = we have invariant surfaces £3 = 1 + ax and £4 = l + ax — ey/2 
with cofactors ax and ax + ey, giving a hnearizing change of coordinates 

{x,Y,z) = {xq\y4e^\zq%). 

Case 17: This case has invariant algebraic surfaces £1 = 1 + ax — ^y = and 
£2 = 1 + kz + ^xz — bkyz = 0, with cofactors = ax + by and Li^ = kz. 
From these we obtain first integrals (f)i = x"^ y £^^ and (f>2 = y z'^ £2^ , and a 
linearizing change of coordinates {X,Y,Z) = {x£^ ,y£i ,z£i£2 )• 

Case 18: In this case the system has an invariant algebraic plane £ = 
(1 + ax + by)'^ + kz{l + |x — by) = with cofactor Li = 2ax — 4by + kz. It is 
easy to obtain two independent first integrals (pi = x'^ y and <t>2 = y z"^ £~'^ ■ 

Case 19: In this case we find two invariant algebraic surfaces = (1 + 
ax + byY — acxz = and £2 = {1 + ax + by)"^ — 2cz(l + ax — by — ^z) = 
with cofactors L^j = 2ax — Aby and Li^ = 2ax — Aby — 2cz. These give the 
first integrals 0i = y£i^ £2 and (f>2 = y z"^ £^^£2^ ■ 

Case 20: The system ([1]) reduces to 

x = x{l + ax + cz), y = y{-2 + ey + fz), z = z{l-fz), 

The first and third equations gives a linearizable node and the change of 
coordinates Y = fiV'^l linearizes the second equation. 

Case 21: The system in this case has three invariant algebraic surfaces 
£1 = I + ax - ^y = 0, £2 = I + ax = and £3 = I - §y + fxy = with 
cofactors Li-^ = ax + ey, L^,^ = ax and = ey. It is easy to find two 

independent first integrals cpi = x'^y£^^ and (f)2 = y £2°' ^3 " ! and a 
linearizing change of coordinates (X, Y, Z) = {x £2 -\y£2£3\z£2^£3 =). 

When a = 0, we replace ^2^" with the exponential factor exp(x), while, 
when e = 0, we replace l^'^ by the exponential factor exp(— y(l — ax)/2). 

□ 
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